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[i]  This  study  examines  the  contribution  of  ocean  dynamics  to  sea  surface  temperature 
(SST)  biases  in  the  eastern  Pacific  cold  tongue  region  in  fifteen  coupled  general  circulation 
models  (CGCMs)  participating  in  the  Intergovernmental  Panel  on  Climate  Change  (IPCC) 

Fourth  Assessment  Report  (AR4).  Twenty  years  (1980-1999)  of  the  twentieth-century 
(20C3m)  climate  simulations  from  each  model  are  analyzed.  An  excessive  and  narrow  SST 
cold  tongue  that  extends  too  far  west  into  the  western  Pacific  in  comparison  to  observations 
is  a  common  bias  in  CGCMs.  This  feature  is  found  in  CGCMs  analyzed  here  and  in 
many  previous  studies.  The  heat  budget  analysis  indicates  that  errors  in  both  net  surface  heat 
flux  and  total  upper  ocean  heat  advection  significantly  contribute  to  the  excessive  cold 
tongue  in  the  equatorial  Pacific.  The  stronger  heat  advection  in  the  models  is  caused  by 
overly  strong  horizontal  heat  advection  associated  with  too  strong  zonal  currents,  and  overly 
strong  vertical  heat  advection  due  to  excessive  upwelling  and  the  vertical  gradient  of 
temperature.  The  Bjerknes  feedback  in  the  coupled  models  is  shown  to  be  weaker  than  in 
observations,  which  may  be  related  to  the  insufficient  response  of  surface  zonal  winds  to 
SST  in  the  models  and  an  erroneous  subsurface  temperature  structure.  A  hypothesis  that 
describes  how  the  cold  tongue  bias  is  possibly  developed  in  the  CGCMs  is  provided  based 
on  the  results  of  our  analysis. 

Citation:  Zheng,  Y.,  J.-L.  Lin,  and  T.  Shinoda  (2012),  The  equatorial  Pacific  cold  tongue  simulated  by  IPCC  AR4  coupled 
GCMs:  Upper  ocean  heat  budget  and  feedback  analysis,  J.  Geophys.  Res.,  117 ,  C05024,  doi:  10. 1029/201 1JC007746. 


1.  Introduction 

[2]  The  equatorial  Pacific  is  observed  to  have  a  minimum 
sea  surface  temperature  (SST)  that  extends  from  the  west 
coasts  of  the  Americas  into  the  central  Pacific.  This  extension 
of  cool  water  is  commonly  referred  to  as  the  cold  tongue 
[ Wyrtki,  1981].  It  is  generally  argued  that  the  Pacific  cold 
tongue  is  maintained  by  horizontal  advection  of  cold  water 
from  the  east  and  by  upwelling  of  cold  water  from  the 
subsurface.  Net  heat  flux  at  the  sea  surface  warms  the  two 
contributions.  The  contributions  of  cool  advection  and  of 
upwelling  as  well  as  the  net  surface  heat  flux  vary  throughout 
the  year  and  determine  the  annual  and  interannual  variations 
of  the  cold  tongue  intensity.  Previous  studies  discussed  the 
relative  roles  of  these  contributing  terms  using  large  regional 
box  models  [Wyrtki,  1981;  Bryden  and  Brady,  1985].  For 
example,  Wyrtki  [1981]  suggested  from  the  heat  budget 
analysis  that  although  horizontal  advection  in  the  South 
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Equatorial  Current  (SEC)  contributes  to  the  cold  tongue,  the 
contribution  from  upwelling  is  much  larger.  He  further  indi¬ 
cated  that  upwelling  is  limited  to  relatively  shallow  depths 
above  the  core  of  the  Equatorial  Undercurrent  (EUC)  and  the 
temperatures  are  around  3°C  less  than  the  temperature  of 
the  water  flowing  out  laterally.  However,  Bryden  and  Brady 
[1985]  found  that  the  zonal  and  meridional  convergences 
contribute  equally  to  the  cold  tongue.  They  also  concluded 
that  although  the  upwelling  is  large,  it  contributes  little 
because  it  occurs  in  conjunction  with  the  eastward  flow  of  the 
Equatorial  Undercurrent,  which  shoals  to  the  east  and  flows 
predominantly  along  isotherms.  Using  data  from  satellite- 
tracked  drifting  buoys  and  VOS/XBT  profiles,  Swenson  and 
Hansen  [1999]  estimated  the  seasonal  cycle  of  the  tropical 
Pacific  Ocean  mixed  layer  heat  budget  in  the  cold  tongue 
region  and  concluded  that  oceanic  heat  advection  and  eddy 
heat  flux  convergence  play  a  significant  role  in  the  evolution 
of  the  Pacific  cold  tongue. 

[3]  Besides  the  contribution  of  oceanic  physical  pro¬ 
cesses  as  discussed  in  the  above  studies,  ocean-atmosphere 
feedback  also  plays  an  important  role  in  the  formation 
and  maintenance  of  the  Pacific  cold  tongue.  Bjerknes  [1969] 
first  hypothesized  a  positive  feedback  of  tropical  ocean- 
atmosphere  interaction  that  can  amplify  a  perturbation  of 
SST  in  the  cold  tongue.  The  trade  winds  generate  the  equa¬ 
torial  upwelling,  effectively  bringing  cold  water  from  the 
subsurface  to  the  surface  to  form  a  cold  tongue  in  the  eastern 
Pacific.  Meanwhile,  the  zonal  atmospheric  pressure  gradient 
produced  by  the  zonal  gradient  of  the  SST  drives  the  Walker 
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circulation,  which  intensifies  the  surface  easterly  winds  over 
the  Pacific  Basin,  further  strengthening  the  cold  tongue. 
In  this  study,  we  will  revisit  this  feedback  in  climate  models. 

[4]  It  is  well  known  that  the  Pacific  cold  tongue  plays  an 
important  role  in  the  formation  of  tropical  cloud  and  precip¬ 
itation  patterns,  in  supplying  nutrients  for  surface  ocean 
biological  productivity,  and  in  carbon  cycling  since  it  is  the 
major  oceanic  source  of  C02  for  the  atmosphere  [Field  et  al, 
1998;  Calvo  et  al. ,  2011].  Large  SST  anomalies  associated 
with  El  Nino-Southern  Oscillation  (ENSO)  are  manifested  as 
variations  about  the  mean  state  of  the  Pacific  cold  tongue, 
which  affect  much  of  the  global  climate  and  weather  patterns. 
Despite  its  importance  in  global  hydrological  and  biogeo- 
chcmical  cycles,  as  well  as  its  influences  on  global  weather 
and  climate  patterns,  the  equatorial  Pacific  cold  tongue 
has  not  been  well  simulated  by  coupled  general  circulation 
models  (CGCMs)  used  for  climate  predictions  and  projec¬ 
tions  [Neelin  et  al ,  1992;  Mechoso  et  al ,  1995;  Delecluse 
et  al ,  1998;  Latif  et  al ,  2001;  Davey  et  al ,  2002;  Meehl 
et  aly  2005;  Wittenberg  et  al ,  2006;  Dai ,  2006;  Lin,  2007]. 
For  example,  in  comparison  to  observations,  an  excessive 
and  overly  narrow  sea  surface  cold  tongue  that  extends  too 
far  west  into  the  western  Pacific  in  the  couple  climate  models 
has  been  a  common  feature  for  almost  two  decades  [Mechoso 
et  al,  1995;  Latif  et  al.y  2001;  Randall  et  aly  2007].  This 
equatorial  Pacific  cold  bias  is  usually  associated  with  the 
double-intcrtropical  convergence  zone  (ITCZ)  problem,  which 
is  characterized  by  excessive  precipitation  off  the  equator  and 
insufficient  precipitation  on  the  equator  in  comparison  with 
observations. 

[5]  The  origin  of  the  Pacific  cold  tongue  bias  in  climate 
models  is  complex.  Several  hypotheses  have  been  proposed 
to  understand  the  possible  causes  of  the  Pacific  cold  tongue 
bias.  Those  possible  causes  include  a  wrong  deep  convection 
feedback  with  the  planetary  boundary  layer  (PBL)  and  the 
SST,  wrong  solar  and  longwave  cloud  forcing,  bad  coastal 
winds  resulting  from  coarse  resolutions  and  continental 
effect,  wrong  ocean  mixing  schemes,  wrong  air-sea  inter¬ 
actions  giving  rise  to  an  improperly  stabilized  PBL  over  the 
cold  tongue,  large  coupling  intervals  between  atmospheric 
GCM  and  oceanic  GCM,  and  no  biology  embedded  to  affect 
penetrative  radiation  in  climate  models  (P.  Chang,  Summary 
of  the  Tropical  Bias  Workshop,  May  28-30,  2003,  GFDL, 
Princeton  Univ.,  Princeton,  N.  J.,  available  at  ftp://grads.iges. 
org/pub/schneider/CTBW05/Previous_mtg/Ping_CLIVAR_ 
summary. pdf)-  Consequently,  there  have  been  several  efforts 
to  reduce  the  SST  bias  in  the  Pacific  cold  tongue  in  CGCMs. 
For  example,  Yu  and  Mechoso  [1999]  reduced  this  cold 
tongue  bias  by  imposing  an  annually  varying  stratus  cloud 
amount  off  the  Peruvian  coast.  The  bias  can  also  be  reduced 
by  refining  horizontal  and  vertical  resolutions  (C.  R.  Mechoso, 
Modeling  the  southeastern  Pacific  climate:  Progress 
and  challenges,  NCEP  EMC  seminar,  2006;  available  at 
http://www.emc.ncep.noaa.gov/seminars/presentations/2006/ 
Mechoso.NCEP.Jan.06.ppt),  improving  the  convection 
scheme  [Frey  et  al.y  1997;  Zhang  and  Wang ,  2006;  Song 
and  Zhangy  2009;  Zhang  and  Song ,  2010;  Bellucci  et  al.y 
20 1 0],  and/or  changing  surface  wind  stress  formulation  [Luo 
et  aly  2005]  in  atmospheric  models.  It  is  also  believed  that 
ocean  atmosphere  feedback  plays  a  critical  role  in  determin¬ 
ing  the  tropical  mean  climate.  The  equatorial  Pacific  cold  bias 
can  be  moderately  reduced  by  reducing  the  coupling  interval 


to  three  hours  in  a  coupled  climate  model  [Misra  et  aly  2008]. 
The  ocean-atmosphere  feedback  in  the  tropical  mean  climate 
system  includes  the  SST  gradient-trade  wind  feedback  or 
Bjerknes  feedback  [e.g.,  Bjerkncsy  1969],  the  SST-surface 
latent  heat  flux  feedback  [e.g.,  Wallacey  1992],  and  the  SST- 
surface  shortwave  flux  feedback  [e.g.,  Ramanathan  and 
CollinSy  1991;  Peters  and  Brethertony  2005].  In  fact,  many 
theories  have  been  developed  for  ocean-atmosphere  feedback 
mechanisms  in  the  coupled  tropical  mean  climate  [Neelin  and 
Dijkstray  1995;  Dijkstra  and  Neelin y  1995,  1999;  Sun  and  Liuy 
1996;  Jin,  1996;  Clement  et  al.y  1996,  2005;  Liu ,  1997;  Liu 
and  Huangy  1997;  Clement  and  Seager,  1999;  van  der  Vaart 
et  aly  2000;  Cal  2003]. 

[6]  Although  these  previous  studies  have  improved  our 
understanding  of  the  mean  climate  problem  including  the 
equatorial  Pacific  cold  tongue  bias,  many  issues  remain 
unresolved.  For  example,  the  contribution  of  ocean  heat 
advection  to  the  equatorial  Pacific  cold  tongue  bias  in  climate 
models  is  still  unknown.  For  the  Bjerknes  feedback,  previous 
studies  [e.g.,  Lin ,  2007]  have  examined  only  the  strength  of 
its  atmospheric  component  (i.e.,  response  of  winds  to  SST 
gradient);  none  have  investigated  the  strength  of  its  oceanic 
component  (i.e.,  response  of  ocean  currents  and  heat  advec¬ 
tion  to  surface  winds)  or  the  strength  of  the  entire  feedback 
loop,  both  of  which  require  the  calculation  of  the  upper  ocean 
heat  budget.  Furthermore,  Lin  [2007]  evaluated  the  response 
of  winds  to  SST  using  the  atmospheric  model  simulations 
forced  by  observed  SST  to  isolate  the  possible  biases  in 
atmospheric  models.  However,  when  the  atmospheric  mod¬ 
els  are  coupled  with  ocean  models,  the  wind  response  to  SST 
may  be  changed  by  any  biased  mean  state  or  by  coupling.  In 
this  study,  we  will  focus  on  the  equatorial  Pacific  cold  tongue 
bias  by  examining  the  fully  coupled  simulations  of  15  IPCC 
AR4  CGCMs  and  directly  calculating  the  upper  ocean  heat 
budget  and  feedback  strength.  The  purposes  are  threefold* 

[7]  1.  Examine  the  structure  of  ocean  circulation  and 
thermodynamics  associated  with  the  equatorial  Pacific  cold 
tongue  bias  in  CGCMs. 

[8]  2.  Examine  whether  ocean  heat  advection  significantly 
contributes  to  local  SST  bias  on  the  equator. 

[9]  3.  Examine  the  strength  of  the  entire  Bjerknes  feedback 
loop  and  its  oceanic  component  in  CGCMs. 

[10]  The  remainder  of  this  paper  is  organized  as  follows. 
The  models  and  validation  data  sets  used  in  this  study  are 
described  in  section  2.  The  structure  of  upper  ocean  circula¬ 
tion  in  the  tropical  Pacific  is  presented  in  section  3.  In  section 
4  the  contribution  of  ocean  heat  advection  to  SST  bias  in  the 
equatorial  Pacific,  which  is  based  on  an  analysis  of  the  upper 
ocean  heat  budget,  is  evaluated.  In  section  5,  the  Bjerknes 
feedback,  including  its  oceanic  component,  is  discussed. 
Finally,  a  summary  is  presented  in  section  6. 

2.  Models  and  Validation  Data  Sets 
2.1.  IPCC  Models 

[11]  The  analysis  is  based  on  20-ycar  (1980-1999)  model 
runs  of  the  Climate  of  the  Twentieth  Century  (20C3M) 
simulations  from  15  coupled  GCMs.  Table  1  shows  the 
model  names  and  acronyms,  their  ocean  model  horizontal 
and  vertical  resolutions,  the  heat  flux  corrections,  and  the 
runs  employed  for  analysis.  For  each  model,  we  use  20  years 
of  monthly  mean  ocean  temperature,  three-dimensional  ocean 


2  of  20 


C05024 


ZHENG  ET  AL.:  THE  PACIFIC  COLD  TONGUE  BIAS  ANALYSIS 


C05024 


Table  1.  List  of  15  IPCC  AR4  Coupled  GCMs  That  Participated  in  This  Study 


Modeling  Groups 

IPCC  ID 

(Label  in  Figures) 

Resolution* 

Heat  Flux 
Correction 

Run 

Number 

Canadian  Centre  for  Climate  Modeling  and  Analysis 

CGCMA3.1-T47  (cgcm-t47) 

192  x  96-L29 

Yes 

I 

Canadian  Centre  for  Climate  Modeling  and  Analysis 

CGCMA3.1-T63  (cgcm-t63) 

256  x  I92-L29 

Yes 

1 

M£t6o  France  /  Centre  National  dc  Rechcrchcs  M&corologiqucs 

CNRM-CM3  (enrm) 

180  x  1 70-L33 

None 

1 

CSIRO  Atmospheric  Research 

CS1RO-Mk3.0  (csiro-mk3.0) 

192  x  I89-L3I 

None 

2 

CS1RO  Atmospheric  Research 

CSIRO-Mk3.5  (csiro-mk3.5) 

192  x  189-L31 

None 

1 

NASA/Goddard  Institute  for  Space  Studies 

GISS-AOM  (giss-aom) 

90  x  60-L31 

None 

1 

NASA/Goddard  Institute  for  Space  Studies 

G1SS-ER  (giss-cr) 

72  x  46-L33 

None 

1 

LASG/lnstitutc  of  Atmospheric  Physics 

FGOALS-gl.O  (iap) 

360  x  170-L33 

None 

1 

Instituto  Nazionalc  di  Gcofisica  c  Vulcanologia 

INGV-SXG  (ingv) 

360  x  I80-L33 

None 

1 

Institute  Pierre  Simon  Laplace 

1PSL-CM4  (ipsl) 

180  x  I70-L3I 

None 

I 

Center  for  Climate  System  Research  (The  University  of  Tokyo), 
National  Institute  for  Environmental  Studies,  and  Frontier 
Research  Center  for  Global  Change 

MIROC3.2  -medres  (miroc-mcdrcs) 

256  x  192-L33 

None 

2 

Max  Planck  Institute  for  Meteorology 

ECHAM5/MPI-OM  (mpi) 

360  x  180-L40 

None 

1 

Meteorological  Research  Institute 

MRJ-CGCM2.3.2  (mri) 

144  x  1II-L23 

Yes 

1 

National  Center  for  Atmospheric  Research 

CCSM3  (ccsm3) 

320  x  395-L40 

None 

1 

Hadley  Centre  for  Climate  Prediction  and  Rcscarch/Mct  Office 

UKMO-HadGEM  1  (hadgcml) 

360  x  216-L40 

None 

1 

“Resolution  is  about  ocean  model  output  which  is  denoted  by  grid  points  in  longitude  x  latitude  and  number  of  vertical  layers. 


currents,  surface  wind  stress,  surface  downward/upward 
shortwavc/longwavc  radiation,  surface  latent  heat  flux,  and 
surface  sensible  heat  flux. 

2.2.  SODA 

[12]  The  Simple  Ocean  Data  Assimilation  (SODA)  meth¬ 
odology,  the  ingested  data,  and  the  error  covariance  structure 
of  both  the  model  and  the  observations  arc  described  by 
Carton  et  al.  [2000a,  2000b],  Carton  and  Giese  [2008],  and 
Zheng  and  Giese  [2009].  In  this  study,  SODA  v2.0.2  is  used. 
SODA  ocean  analysis  is  suitable  for  evaluating  the  simula¬ 
tions  of  CGCMs  without  assimilating  data  [Zheng  et  al ., 
2011].  The  ocean  model  is  based  on  the  Los  Alamos  imple¬ 
mentation  of  the  Parallel  Ocean  Program  (POP)  [Smith  et  al., 
1992].  The  model  resolution  is  on  average  0.4°  (Ion)  x 
0.25°  (lat)  with  40  levels  in  the  vertical.  The  model  is  forced 
with  the  European  Centre  for  Medium-Range  Weather 
Forecasts  (ECMWF)  ERA-40  daily  atmospheric  reanalysis 
winds  [Simmons  and  Gibson ,  2002]  for  the  44-year  period 
from  1958  to  2001. 

[13]  The  model  is  constrained  by  abundant  temperature 
and  salinity  observations  using  a  sequential  assimilation 
algorithm,  which  is  described  by  Carton  et  al.  [2000a, 
2000b]  and  Carton  and  Giese  [2008].  Surface  heat  fluxes  are 
computed  from  bulk  formulae  [Smith  et  al .,  1992],  using 
atmospheric  variables  from  the  NCEP/NCAR  reanalysis 
[Kalnay  et  al. ,  1996].  The  NCEP/NCAR  reanalysis  infor¬ 
mation  is  used  for  the  bulk  formulae  instead  of  the  ERA-40 
variables  throughout  the  experiment  to  give  continuity  of 
surface  forcing  during  periods  for  which  the  ERA-40  winds 
are  not  available.  However,  the  details  of  the  surface  heat  flux 
boundary  condition  are  relatively  unimportant  in  influencing 
the  solution,  since  near-surface  temperature  observations  are 
used  to  update  the  mixed  layer  temperature.  Vertical  diffu¬ 
sion  of  momentum,  heat,  and  salt  is  based  on  a  nonlocal 
K-profilc  parameterization  (KPP)  [Large  et  al. ,  1994]  and 
horizontal  diffusion  for  subgrid-scalc  processes  is  based  on  a 
biharmonic  mixing  scheme. 

[14]  Averages  of  model  output  variables  (temperature, 
salinity,  and  velocity)  arc  saved  at  five-day  intervals.  These 
average  fields  are  remapped  onto  a  uniform  global  0.5°  x  0.5° 
horizontal  grid  using  the  horizontal  grid  spherical  coordinate 


remapping  and  interpolation  package  with  second-order  con¬ 
servative  remapping  [Jones,  1999]. 

2.3.  Heat  Flux  Data  Set 

[15]  In  this  study,  monthly  mean  surface  heat  fluxes 
from  the  Objectively  Analyzed  air-sca  Heat  Fluxes  Project 
(OAFlux)  [Yu  and  Weller ,  2007;  Yu  et  al .,  2008]  at  the 
Woods  Hole  Oceanographic  Institution  (WHOl)  arc  used  for 
evaluating  the  heat  fluxes  in  CGCMs  because  they  are  the 
latest,  and  perhaps  the  best  validated  data  sets.  Daily  fluxes 
are  computed  from  the  optimally  estimated  variables  using 
the  state-of-the-art  Tropical  Ocean  Global  Atmosphere 
Coupled  Ocean-Atmosphere  Response  Experiment  (TOGA 
CO  ARE)  bulk  flux  algorithm  3.0  [Fairall  et  al. ,  2003]. 
Surface  shortwave  and  longwave  radiation  of  OAFlux  is 
derived  from  the  International  Satellite  Cloud  Climatology 
Project  (ISCCP)-FD  estimates  [Zhang  et  al.,  2004],  which 
are  available  from  1  July  1983.  These  flux  estimates  are 
compared  to  107  (105  buoys  and  2  ships)  in  situ  flux  time 
scries,  and  it  is  found  that  they  arc  relatively  unbiased  and 
have  the  smallest  mean  error  compared  to  other  data  sets 
[Yu  and  Weller,  2007;  Yu  et  al.,  2008]. 

3.  Structure  of  Upper  Ocean  Circulation 
and  Thermodynamics 

[16]  First  we  look  at  the  structure  of  upper-ocean  circula¬ 
tion  and  thermodynamics  associated  with  the  equatorial  cold 
tongue  bias  in  the  tropical  Pacific.  To  qualitatively  depict 
the  bias  in  the  position  of  the  equatorial  Pacific  cold  tongue, 
the  observed  and  simulated  ocean  mean  state  of  horizontal 
distribution,  rather  than  the  difference  between  models  and 
observations,  is  shown.  The  cold  tongue  bias  in  magnitude  is 
also  very  important  to  be  quantified,  as  will  be  discussed  in  a 
later  section.  Figure  1  shows  the  spatial  distribution  of  time- 
mean  SST  and  surface  currents  in  the  tropical  Pacific  from 
SODA  ocean  analysis  and  the  15  IPCC  AR4  CGCMs.  Errors 
in  ocean  surface  circulation  and  SST  spatial  patterns  in 
models  are  prominent.  For  example,  compared  to  observa¬ 
tions,  the  cold  tongue  in  most  models  is  too  strong  and 
narrow,  and  it  extends  too  far  west  into  the  western  Pacific. 
Also  there  are  strong  westward  currents  in  the  central  and 
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Figure  1.  Spatial  distribution  of  time-mean  sea  surface  temperature  (SST;  shading  and  contours  in  °C) 
and  surface  ocean  currents  (u,v;  vectors  in  m  s“l)  from  SODA  ocean  analysis  and  15  IPCC  AR4  CGCMs 
in  the  tropical  Pacific  (120°E-70°W,  20°S-20°N).  The  time-mean  values  arc  computed  by  averaging 
monthly  data  over  January  1980-Dcccmbcr  1999.  Contour  interval  for  SST  is  1°C. 


cast  equatorial  Pacific.  The  cold  SST  biases  in  the  equatorial 
central  and  eastern  Pacific  may  be  associated  with  stronger- 
than-observed  westward  surface  currents  and  their  associated 
heat  advection  along  the  equator,  in  addition  to  coupled 
processes  such  as  cloud-radiation,  wind-ovaporation  feed¬ 
back  as  discussed  in  detail  by  Lin  [2007]. 

[n]  Figure  2  shows  the  spatial  distribution  of  annual  mean 
temperature  and  ocean  currents  averaged  in  the  upper  50  m  in 
the  tropical  Pacific  Ocean.  The  spatial  pattern  in  temperature 
for  each  model  is  similar  to  that  for  SST  in  Figure  1  but  with  a 
slightly  smaller  magnitude  due  to  depth  averaging;  thus,  the 
signature  of  the  cold  tongue  bias  found  in  the  ocean  surface  is 
also  present  in  the  upper  ocean.  In  the  equatorial  Pacific, 
most  models  produce  stronger  westward  currents  than 
observed  that  transport  water  from  the  cold  tongue  region  to 
the  western  warm  pool  region.  Therefore,  they  may  contrib¬ 
ute  to  cold  SST  bias  through  overly  strong  cold  advection 
in  addition  to  equatorial  upwelling.  Thus,  the  far  westward 
extension  of  the  cold  tongue  in  most  models  is  associated 


with  overly  strong  westward  surface  and  subsurface  currents 
that  efficiently  bring  cold  water  far  westward.  Large  biases 
in  upper  ocean  circulation  in  giss-aom  and  giss-er  arc  one  of 
the  main  causes  of  large  warm  SST  biases  in  the  entire 
equatorial  Pacific. 

[is]  Equatorial  upwelling  is  one  of  the  major  processes 
affecting  equatorial  SST.  Similar  to  Figure  1,  ocean  mean 
state  on  a  vertical  plane  is  shown  to  better  qualitatively 
describe  the  biases  of  actual  thermocline  position  and  actual 
vertical  circulation  in  CGCMs  compared  to  observations. 
Figure  3  illustrates  the  vertical  cross  section  of  annual  mean 
temperature  and  ocean  circulation  averaged  between  2°S 
and  2°N.  The  black  curve  is  the  20°C  isotherm,  which  is  in 
the  middle  of  the  equatorial  main  thermocline.  The  eastward 
EUC  in  the  central  and  east  equatorial  Pacific  is  too  weak  in 
most  CGCMs.  Models  also  have  problems  simulating  the 
vertical  profile  of  upper  ocean  temperature.  For  example,  the 
deepest  equatorial  thermocline  in  some  models  (e.g.,  egem- 
t47 ,  cgcm-t63 ,  csiro-mk3.0,  csiro-mk3.5 ,  miroc-medres ,  mpi) 
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Figure  2.  As  in  Figure  1,  but  for  those  averaged  in  the  upper-50-m  layer. 


is  not  at  the  far  end  of  the  western  Pacific  (about  130CE  in 
SODA),  but  between  140°E  and  160°E.  Similarly,  the  shal¬ 
lowest  equatorial  thermocline  in  some  models  is  not  at  the 
east  coastline  (around  80°  W  in  SODA)  but  at  around  100°W 
(e.g.,  csiro-mk3.0 ,  csiro-mk3.5 ,  iap,  ingv,  ipsl,  miroc-medres, 
mpi,  mri,  ccsm3)  or  around  120°W  (e.g.,  hadgeml).  Since  the 
thermocline  shoals  in  the  eastern  Pacific,  errors  in  vertical 
velocity  can  produce  significant  SST  biases.  The  upwelling  is 
too  strong  between  140°E  and  160°W  in  models  such  as 
cnrm,  iap ,  ingv,  mpi  and  hadgeml,  which  partly  explains  the 
far  westward  extension  of  the  cold  tongue.  We  will  examine 
upper  ocean  processes  that  contribute  to  errors  in  the  tem¬ 
perature  profile  and  the  circulation  simulated  in  these  models, 
particularly  in  the  central  and  eastern  Pacific  Ocean.  Since 
temperature  in  the  upper  50  m  is  a  good  proxy  to  time-mean 
SST,  in  this  study  we  focus  our  analysis  on  the  ocean  pro¬ 
cesses  in  this  layer. 

4.  Upper  Ocean  Heat  Budget  Analysis 

[19]  In  this  section,  we  analyze  the  upper  ocean  heat  bud¬ 
get  to  understand  oceanic  processes  that  could  contribute  to 


the  SST  bias  in  IPCC  CGCMs.  The  method  for  calculating 
the  upper  ocean  heat  budget  follows  that  of  Zheng  et  al 
[2010,  2011].  Note  that  CGCMs  are  not  eddy  resolved; 
thus,  the  eddy  heat  flux  divergence  term  is  not  calculated. 

4,1,  Net  Surface  Heat  Flux 

[20]  Figure  4  shows  the  horizontal  distribution  of  net  sur¬ 
face  heat  flux  in  the  entire  Pacific  from  OAFIux  and  CGCMs 
that  were  averaged  over  the  period  July  1983-December 
1999  (positive  values  denote  warming  of  the  ocean).  Note 
that  ingv  is  not  included  here  because  surface  shortwave  and 
longwave  fluxes  are  not  available  for  this  model.  Also,  note 
that  radiation  from  ISCCP-FD  (and  thus  net  surface  fluxes 
from  OAFIux)  is  available  only  from  1  July  1983  for  the 
comparison.  Most  models  produce  weaker  (i.e.,  less  heat 
input  to  the  ocean)  net  surface  heat  fluxes  compared  to 
OAFIux  in  the  entire  tropical  Pacific.  The  belt  of  large  sur¬ 
face  heat  flux  in  the  central  and  eastern  Pacific  is  too  narrow 
in  most  models,  compared  with  the  observation,  except  for 
the  three  models  ( cgcm-t47 ,  cgcm-t63,  and  mri)  for  which 
heat  flux  correction  was  applied.  A  weak  cooling  resulting 
from  surface  heat  flux  is  present  over  some  regions  off  the 
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Figure  3.  Annual  mean  vertical  ocean  circulation  (vectors  inms"1)  and  temperature  (shading  contours 
in  °C)  in  the  zonal-vertical  plane  averaged  along  2°S-2°N  from  monthly  data  in  SODA  ocean  analysis 
and  15  IPCC  AR4  CGCMs.  20°C  isotherm  is  denoted  by  the  black  curve.  The  period  for  averaging  is 
January  1 980-Dcccmbcr  1999. 


equator  in  all  models,  whereas  in  OAFlux  there  is  warming 
over  the  entire  tropic  region.  These  errors  in  the  net  surface 
heat  flux  contribute  to  the  SST  cold  biases  along  the  equator. 

4.2.  Heat  Advection 

[21]  Total  heat  advection  integrated  from  the  surface  to  a 
depth  zq  is  defined  as 

*0  r  0  /  3  «r»  j>r\ 

V  ’  VTdz  "  -  i  pC\u^  +  +  w 

where  p  is  the  density  of  seawater;  Cp  is  the  specific  heat 
capacity  of  seawater  at  constant  pressure;  V  is  the  3D 
velocity  vector,  which  is  u\  +  vj  +  uk;  and  w,  v,  w  are 
their  components  in  zonal,  meridional,  and  vertical  direc¬ 
tion,  respectively.  The  3D  gradient  vector,  V,  is  defined  as 
d!dx\  +  d!dy\  +  dldzk.  Total  heat  advection  in  the  upper  50  m 


from  the  15  IPCC  AR4  CGCMs  is  calculated  for  the  period 
January  1 980-December  1999,  and  is  then  compared  with 
the  total  heat  advection  from  SODA. 

[22]  Figure  5  shows  the  spatial  distribution  of  temporal 
mean  total  heat  advection  from  the  15  IPCC  AR4  CGCMs 
and  SODA.  The  total  heat  advection  in  the  upper  50  m 
contributes  to  cooling  in  the  equatorial  Pacific  Ocean  in  the 
models  and  in  SODA.  In  11  of  15  models  (all  except  egem- 
t47 ,  cgcm-t63y  giss-aom ,  and  giss-er),  the  upper-50-m  total 
heat  advection  contributes  more  cooling  to  the  equatorial 
Pacific  Ocean  than  SODA  does.  Interestingly,  in  these  spe¬ 
cific  models,  an  excessive  and  narrow  total  heat  advection 
area,  extending  far  west  into  the  western  Pacific,  is  similar  to 
the  shape  of  the  cold  tongue  simulated  in  the  models.  The 
excessive  cooling  from  total  heat  advection  and  the  insuffi¬ 
cient  warming  from  net  surface  heat  flux  in  most  models  act 
together  to  produce  SST  cold  biases  in  the  equatorial  Pacific. 
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Figure  4.  Temporal  mean  net  surface  heat  flux  (W  m  2)  in  the  tropical  Pacific  Ocean  from  OAFlux  and 
14  IP CC  AR4  CGCMs.  The  period  for  averaging  is  July  1983-December  1999. 
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Figure  5.  Annual  mean  total  heat  advection  (W  m  2)  in  the  upper-50-m  layer  in  the  tropical  Pacific  Ocean 
from  SODA  ocean  analysis  and  15  IPCC  AR4  CGCMs.  The  period  for  averaging  is  January  1980-December 
1999. 
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Figure  6.  Annual  mean  (a)  SST  (°C),  (b)  net  surface  heat  flux  (W  m“2),  and  (c)  upper-50-m  total  heat 
advection  (W  m“2)  averaged  along  2°S-2°N  from  observed  monthly  data  and  15  IPCC  AR4  CGCMs. 
The  data  are  first  converted  to  10-degree  bins  in  longitude  before  meridional  averaging.  The  observed  data 
for  SST,  net  surface  heat  flux,  and  zonal  wind  stress  are  from  SODA  ocean  analysis,  OAFlux,  and  ERA-40, 
respectively.  The  period  for  averaging  is  January  1 980-December  1999,  except  for  OAFlux,  which  is  aver¬ 
aged  over  the  period  January  1 984-Dccember  1999. 


Therefore,  errors  in  heat  advection  contribute  to  the  SST  cold 
biases  in  the  equatorial  Pacific.  It  should  be  noted  that  the 
large  error  in  total  heat  advection  from  giss-aom  and  giss-er 
is  primarily  due  to  the  poor  simulation  of  upper  ocean  cir¬ 
culation  (Figures  1  and  2). 

[23]  Figure  6  displays  longitudinal  variations  of  SST,  net 
surface  heat  flux,  and  total  heat  advection  that  are  averaged 
over  the  region  of  2°S-2°N  for  January  1 980-December 
1999,  except  for  OAFlux,  which  is  averaged  over  the  period 
from  January  1984  to  December  1999.  Note  that  the  spatial 
pattern  and  magnitude  of  annual  mean  net  surface  flux  in 
models  averaged  over  the  period  January  1 984- -December 
1999  are  very  similar  to  those  in  Figure  6.  Figure  6  clearly 
indicates  the  quantitative  biases  in  cold  tongue  SST,  net 


surface  heat  flux,  and  total  heat  advection  in  climate  models. 
The  unrealistic  westward  extension  of  the  simulated  cold 
tongue  is  clearly  seen  in  Figure  6a.  Most  models  generate 
cold  SST  biases  west  of  120°W  in  the  equatorial  Pacific, 
which  is  associated  with  insufficient  warming  from  net  sur¬ 
face  heat  flux  (Figure  6b)  and  excessive  cooling  from  total 
heat  advection  west  of  120°  W  (Figure  6c). 

4.3.  Relative  Roles  of  Horizontal  and  Vertical  Heat 
Advection 

[24]  Figure  7  shows  the  horizontal  distribution  of  annual 
mean  horizontal  heat  advection  in  the  upper  50  m  of  the 
tropical  Pacific  Ocean.  The  large  horizontal  heat  advection  is 
mostly  found  between  5°S  and  5°N.  In  the  equatorial  Pacific 
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Figure  7.  As  in  Figure  5,  but  for  horizontal  heat  advection. 


there  is  a  dipole-like  structure  in  the  horizontal  heat  advec¬ 
tion,  and  the  weakest  values  are  along  the  equator.  Note  that 
horizontal  heat  advection  causes  warming  in  the  east  equa¬ 
torial  Pacific  adjacent  to  the  coast  (e.g.,  cast  of  100°W). 

[25]  Figure  8  shows  the  horizontal  distribution  of  the 
uppcr-50-m  vertical  heat  advection.  The  large  vertical  heat 
advection  due  to  equatorial  upwelling  is  mostly  found  in  the 
eastern  Pacific  in  observations  and  models.  However,  the 
discrepancies  between  models  and  observations  are  promi¬ 
nent.  For  example,  the  observed  cooling  due  to  strong 
equatorial  upwelling  is  primarily  found  cast  of  140°  W,  while 
11  of  15  models  (all  except  cgcm-t47 ,  cgcm-t63 ,  giss-aom , 
and  giss-er)  have  the  area  of  cooling  extending  too  far  west 
relative  to  observations. 

[26]  To  examine  the  relative  roles  of  vertical  and  horizontal 
heat  advection  in  the  equatorial  Pacific  Ocean,  horizontal  and 
vertical  heat  advection  in  the  CGCMs  is  averaged  over  the 
region  2°S-2°N  and  compared  with  observations  (Figure  9). 
Warming  (up  to  75  W  m~ 2)  due  to  horizontal  heat  advection 
in  the  eastern  equatorial  Pacific  adjacent  to  the  coast  (east 
of  120°W)  is  largely  compensated  for  by  cooling  (up  to 
100  W  m'  2)  due  to  vertical  heat  advection  in  models,  and 
thus,  the  overall  cooling  of  up  to  —25  W  m“2  is  generated  in 


this  region  (Figure  5).  However,  cooling  from  vertical  heat 
advection  in  the  central  and  eastern  equatorial  Pacific  (west 
of  140°W)  is  generally  weaker  than  cooling  from  horizontal 
heat  advection.  Note  that  most  models  have  weaker  vertical 
and  horizontal  heat  advection  near  the  coast  relative  to 
SODA.  Strong  cooling  due  to  horizontal  heat  advection  is 
found  in  the  eastern  and  western  Pacific  whereas  the  cooling 
due  to  vertical  heat  advection  decreases  rapidly  from  east  to 
west.  Therefore,  the  errors  of  both  horizontal  and  vertical 
heat  advection  may  contribute  to  the  unrealistic  cold  tongue 
simulation.  For  example,  the  combined  errors  of  horizontal 
and  vertical  heat  advection  (Figure  6c)  cause  an  excessive 
cooling  up  to  —75  W  m“2  in  the  central  and  eastern  Pacific 
Ocean. 

4.4.  Errors  in  Vertical  and  Horizontal  Heat  Advection 

[27]  Errors  in  vertical  heat  advection  can  be  caused  by 
errors  in  vertical  velocity,  in  vertical  temperature  gradient,  or 
in  both.  Figure  10  shows  the  horizontal  distribution  of  annual 
mean  vertical  velocity  (upward  positive)  averaged  in  the 
upper  50  m  in  SODA  and  in  the  models.  All  models  except 
cgcm-t47 ,  cgcm-t63 ,  giss-aom ,  and  giss-er  generate  a  narrow 
latitudinal  band  (2°S-2°N)  of  strong  upward  velocity.  Nine 
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Figure  8.  As  in  Figure  5,  but  for  vertical  heat  advection. 
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Figure  9.  Annual  mean  (a)  vertical  and  (b)  horizontal  heat  advection  (W  m-2)  in  the  upper  50  m  averaged 
along  2°S-2°N  from  monthly  data  in  SODA  ocean  analysis  and  15  IPCC  AR4  CGCMs.  The  data  are  first 
converted  to  10-degree  bins  in  longitude  before  meridional  averaging.  The  period  for  averaging  is  January 
1980-December  1999. 
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Figure  10.  As  in  Figure  5,  but  for  vertical  velocity. 


of  the  above  1 1  models  (i.e.,  1 1  models  are  all  except  egem- 
t47y  cgcm-t63y  giss-aom ,  and  giss-er)  generate  stronger-than- 
observed  upward  vertical  velocity  and  the  other  two  models 
(iap,  and  ccsm3)  produce  a  slight  weaker  upward  vertical 
velocity  than  observations  either  in  the  eastern  equatorial 
Pacific  (i.e.,  iap)  or  in  the  central  equatorial  Pacific  (i.c., 
ccsm3).  Thus,  overly  strong  upward  velocity  in  most  models 
can  contribute  to  the  SST  cold  biases  of  the  equatorial 
Pacific.  Figure  1 1  shows  the  horizontal  distribution  of  the 
annual  mean  vertical  temperature  gradient  (upward  positive) 
averaged  in  the  upper  50  m  in  SODA  and  in  the  models.  The 
large  vertical  temperature  gradient  is  mostly  confined  to  the 
eastern  Pacific  in  SODA  and  in  most  models  (except  giss- 
aom  and  giss-er).  In  the  equatorial  Pacific,  the  large  vertical 
temperature  gradient  in  six  models  ( csiro-mk3.0y  csiro- 
mk3.5y  ingy ,  ipsl,  miroc-medres ,  and  mpf)  is  generally  located 
cast  of  120°W,  similar  to  the  location  in  SODA.  Five  models 
(cgcm-t47y  cgcm-t63,  iap,  mri,  hadgeml)  have  a  large  verti¬ 
cal  temperature  gradient  extending  west  of  120°W  along  the 
equator.  However,  such  a  large  vertical  temperature  gradient 
and  its  westward  extension  in  cgcm-t47  and  cgcm-t63  do  not 
produce  strong  cold  vertical  heat  advection  because  the 


vertical  velocity  is  very  small  (Figure  10).  Models  such  as 
giss-aom  and  giss-er  produce  weak  vertical  heat  advection 
owing  to  small  vertical  velocity  and  small  vertical  gradient 
of  temperature. 

[28]  Horizontal  heat  advection  near  the  equator  is  mainly 
produced  by  its  zonal  component  (i.e.,  -udT/dx)  since  zonal 
circulation  dominates  along  the  equator.  Figure  12  shows  the 
zonal  profiles  of  annual  mean  zonal  velocity  (m),  vertical 
velocity  (w),  and  zonal  and  vertical  gradient  of  temperature 
averaged  over  the  region  2°S-2°N  and  in  the  0-5  0-m  layer, 
respectively.  It  is  clearly  seen  that  all  models  except  giss-er 
have  westward  currents  that  are  too  strong  in  the  cold  tongue 
region  west  of  120°W  along  with  an  overly  weak  zonal 
gradient  of  temperature.  This  indicates  that  the  excessive 
horizontal  cold  advection  in  the  cold  tongue  region  west  of 
120°W  (Figure  9)  arises  largely  from  the  errors  of  zonal 
currents,  not  from  the  zonal  temperature  gradient.  Eleven 
of  15  models  (all  except  cgcm-t47y  cgcm-t63,  giss-aomy  and 
giss-er)  have  overly  strong  vertical  velocity  as  well  as  an 
overly  large  vertical  temperature  gradient  in  the  cold  tongue 
extension  region.  Therefore,  both  stronger-than-observed 
vertical  velocity  and  the  vertical  gradient  of  temperature  play 
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Figure  11.  As  in  Figure  5,  but  for  the  vertical  gradient  of  temperature. 


an  essential  role  in  producing  excessive  cooling  from  ver¬ 
tical  heat  advection,  which  causes  the  unrealistic  westward 
extension  of  the  cold  tongue. 

5.  Strength  of  the  Bjerknes  Feedback 
in  the  Equatorial  Pacific  Ocean 

[29]  To  discuss  the  Bjerknes  feedback,  we  use  twenty 
years  of  monthly  time  series  which  contain  the  variability  at 
all  time-scales  from  two  months  to  twenty  years  including 
significant  contribution  from  ENSO-related  interannual  var¬ 
iability.  SST  anomaly  can  affect  oceanic  heat  advection 
through  an  atmospheric  feedback  loop.  In  particular,  SST 
anomaly  influences  trade  winds,  which  drive  ocean  currents 
and  upwelling  (downwelling),  thus  affecting  horizontal  and 
vertical  heat  advection.  Heat  advection  in  turn  affects  SST 
directly.  If  there  is  a  linear  relationship  between  SST  and  heat 
advection,  the  SST  anomalies  can  exponentially  grow  or 
decay  depending  on  the  sign  of  the  regression  coefficient.  For 
example,  if  cold  (warm)  advection  occurs  more  often  during 
the  period  of  negative  (positive)  SST  anomalies,  the  cold 
(warm)  SST  anomalies  arc  further  enhanced.  Therefore,  the 
zero-lag  linear  regression  coefficients  of  heat  advection 


versus  SST  can  be  a  measure  of  the  strength  of  feedback 
in  the  coupled  system  in  the  equatorial  Pacific  [e.g.,  Lin, 
2007].  Note  that  the  regression  between  oceanic  heat 
advection  and  SST  represents  the  strength  of  the  entire  loop 
of  the  Bjerknes  feedback,  including  the  atmospheric  response 
to  SST  (not  only  the  oceanic  or  atmospheric  component  of 
the  feedback  loop). 

[30]  Here  we  calculate  the  linear  regression  of  ocean 
heat  advection  versus  temperature  averaged  within  0-50  m 
(T050),  which  is  a  good  proxy  for  mixed  layer  temperature 
(or  SST)  in  this  region.  Figure  1 3  shows  the  linear  regression 
of  2°S-2°N  averaged  of  20-year  monthly  mean  data  for  the 
upper-50-m  vertical  heat  advection  versus  T050  (Figure  13a), 
horizontal  heat  advection  versus  T050  (Figure  13b),  and  total 
heat  advection  versus  T050  (Figure  13c).  We  use  the  same 
vertical  axis  for  the  three  panels  to  compare  the  contribution 
of  vertical  advection  and  horizontal  advection  to  the  total 
advection.  The  analysis  of  observations  indicates  that  a 
positive  feedback  exists  between  vertical  heat  advection 
and  T050  in  the  central  and  eastern  equatorial  Pacific 
Ocean  (Figure  13a),  between  horizontal  heat  advection  and 
T050  (Figure  13b),  and  between  total  heat  advection  and 
T050  (Figure  13c).  This  result  can  be  interpreted  as  follows. 
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Figure  1 2.  Annual  mean  (a)  zonal  velocity  (in  ms  *),  (b)  vertical  velocity  (in  1 0  5  m  s  l),  (c)  horizontal 
temperature  gradient  (in  degK  per  103  km),  and  (d)  vertical  temperature  gradient  (degK  m-1)  in  the  upper 
50  m  averaged  along  2°S-2°N  from  monthly  data  in  SODA  ocean  analysis  and  15  IPCC  AR4  CGCMs. 
The  data  are  first  converted  to  10-degree  bins  in  longitude  before  meridional  averaging.  The  period  for 
averaging  is  January  1 980-December  1999. 


An  increase  in  SST  in  the  eastern  equatorial  Pacific  causes  a 
weaker  zonal  SST  gradient  between  the  cast  and  west  that 
results  in  a  weaker  zonal  wind  stress.  The  weaker  zonal  wind 
stress  generates  a  weaker  equatorial  upwelling  and  zonal 
current,  and  thus  decreases  the  cooling  caused  by  horizontal 
and  vertical  heat  adveetion,  leading  to  a  further  increase  in 
the  SST.  Because  of  this  (Bjerknes)  feedback  loop,  SST 
anomalies  are  enhanced  in  the  coupled  system. 

[31]  An  interesting  result  from  the  analysis  of  observations 
is  that  the  horizontal  heat  adveetion,  which  is  dominated  by 
its  zonal  component,  contributes  more  than  the  vertical 
adveetion  does  to  the  total  adveetion.  This  can  be  understood 
from  the  thermal  structure  of  the  upper  ocean  in  the  central 
and  eastern  equatorial  Pacific  (Figure  3).  In  this  region,  the 
upper  ocean  isotherms  are  tilted  with  longitude:  the  20°C 
isotherm  is  about  1 60  m  deep  at  the  dateline  and  about  50  m 
deep  at  90°W.  The  tilt  of  the  isotherm  means  a  significant 
zonal  temperature  gradient  and  thus  a  substantial  zonal  heat 
adveetion  by  the  zonal  current.  The  traditional  definition  of 
Bjerknes  feedback  considers  only  vertical  heat  adveetion,  but 
Figure  13  shows  that  zonal  heat  adveetion  plays  a  more 


important  role.  Therefore,  the  total  heat  adveetion  should  be 
considered  in  discussions  of  Bjerknes  feedback. 

[32]  Almost  all  models  produce  a  positive  feedback 
between  total  heat  adveetion  and  T050  (Figure  13c),  i.e., 
a  positive  Bjerknes  feedback.  However,  the  models  tend  to 
underestimate  the  strength  of  Bjerknes  feedback  over  the 
eastern  equatorial  Pacific  between  150°W  and  1 10°W,  which 
is  a  key  region  for  the  cold  tongue  and  ENSO.  Biases  in  both 
the  vertical  adveetion  (Figure  13a)  and  horizontal  adveetion 
(Figure  13b)  contribute  to  this  bias.  It  should  be  noted  that  the 
results  are  similar  if  SST,  instead  of  T050,  is  used  to  calculate 
the  regression  (not  shown).  The  overly  weak  Bjerknes  feed¬ 
back  in  CGCMs  is  somewhat  surprising,  as  it  seems  incon¬ 
sistent  with  the  excessive  cold  tongue  in  those  models.  We 
will  discuss  this  issue  further  in  section  6. 

[33]  The  overly  weak  Bjerknes  feedback  in  CGCMs 
could  be  caused  by  many  different  factors  that  are  difficult 
to  isolate  in  the  coupled  system.  Here  we  only  look  at  some 
key  physical  relationships  in  the  models  that  we  hope  we  can 
provide  some  insights  into  the  cause  of  the  problem.  For 
example,  the  regression  between  subsurface  temperature  and 
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Figure  13.  Linear  regression  of  2°S-2°N  averaged  of  20-year  monthly  mean  data  for  (a)  vertical  heat 
advection  (0-50  m)  vs  T050,  (b)  horizontal  heat  advection  (0-50  m)  vs  T050,  (c)  total  heat  advection 
(0-50  m)  vs  T050  (W  m'  2  K_1)  from  monthly  data  in  SODA  ocean  analysis  and  15  IPCC  AR4  CGCMs. 
The  data  are  first  converted  to  10-degree  bins  in  longitude  before  meridional  averaging.  The  period  for  anal¬ 
ysis  is  January  1 980-December  1 999. 
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Figure  14.  Linear  regression  of  2°S-2°N  averaged  of  20-year  monthly  mean  temperature  at  140°W  vs 
SST  from  monthly  data  in  SODA  ocean  analysis,  TAO  buoy  data,  and  15  IPCC  AR4  CGCMs.  The  data 
are  first  converted  to  10-degrce  bins  in  longitude  before  meridional  averaging.  The  period  for  analysis  is 
January  1980-Dcccmbcr  1999.  TAO  buoy  data  from  January  1 985~Deccmbcr  2010  are  used. 


SST  at  14(TW  (Figure  14)  shows  that  subsurface  warming 
(cooling)  is  generally  larger  than  surface  warming  (cooling) 
in  the  observations,  whereas  subsurface  warming  (cooling)  is 
smaller  than  surface  warming  (cooling)  in  most  models.  The 
above  difference  in  subsurface  temperature  structure  between 
models  and  observations  influences  the  relationship  between 
SST,  dT/dz ,  and  vertical  heat  advection,  possibly  giving  rise 
to  an  unrealistic  Bjerknes  feedback.  For  example,  when  the 
same  cold  SST  anomalies  are  found  in  models  and  obser¬ 
vations,  the  cold  anomalies  of  subsurface  temperature  tend  to 
be  smaller  in  models  than  in  the  observations,  leading  to  an 
unrealistic  vertical  temperature  gradient  with  a  sign  opposite 
to  that  in  the  observations.  This  incorrect  vertical  tempera¬ 
ture  gradient  anomaly,  when  combined  with  the  time-mean 
upwelling,  contributes  to  the  weaker  positive  or  even  nega¬ 
tive  vertical  heat  advection  (as  seen  in  Figure  1 3a)  and 
therefore  weakens  the  Bjerknes  feedback.  Figure  15  shows 
the  linear  regression  of  surface  zonal  wind  stress  versus  T050 
(Figure  15a)  and  depth-averaged  vertical  velocity  versus 
surface  zonal  wind  stress  (Figure  15b).  The  analysis  reveals 
that  the  zonal  wind  stress  responds  insufficiently  to  the 
mixed  layer  temperature  (and  thus  SST)  anomaly  in  most 
models  (Figure  15a).  This,  in  combination  with  the  overly 
small  time-mean  zonal  temperature  gradient  in  most  mod¬ 
els  (Figure  12c),  tends  to  produce  insufficient  zonal  heat 
advection  (a  dominant  component  of  horizontal  heat  advec¬ 
tion)  anomaly  associated  with  SST  anomaly  (Figure  13b). 
Therefore,  it  may  contribute  to  the  overly  weak  Bjerknes 
feedback  in  the  models.  As  shown  in  Figure  18a  of  Lin9 s 
[2007]  study,  the  zonal  wind  stress  responds  excessively 
to  the  SST  in  the  atmospheric  models  when  forced  by  the 
observed  SST  (without  coupling  to  an  ocean  model).  The 


weak  response  of  surface  winds  to  T050  (and  SST)  in  the 
coupled  models  suggests  that  coupling  to  an  ocean  model 
and  the  resultant  overly  cold  equilibrium  SST  along  the 
equator  (the  excessive  cold  tongue)  significantly  reduces  the 
excessive  response  of  zonal  winds  to  SST  anomaly. 

[34]  Figure  15b  shows  that  most  models  produce  exces¬ 
sive  upwelling  (downwelling)  anomalies  (i.e.,  more  negative 
values)  in  the  eastern  Pacific  in  response  to  enhanced 
(reduced)  surface  easterly  wind  stress.  This,  in  combination 
with  the  overly  large  time-mean  vertical  temperature  gradient 
in  most  models  (Figure  12d),  tends  to  produce  excessive 
vertical  heat  advection.  Therefore,  the  insufficient  response 
of  vertical  heat  advection  to  mixed  layer  temperature  or  SST 
in  most  of  the  models  shown  in  Figure  13a  is  not  caused  by 
problems  in  the  ocean  dynamical  upwelling  processes  (i.e., 
response  of  ocean  vertical  flow  to  zonal  wind  stress),  but 
more  likely  by  other  processes  such  as  the  insufficient  sur¬ 
face  wind  response  to  SST  anomaly  (Figure  15a)  and  the 
erroneous  subsurface  temperature  structure  (Figure  14). 

6.  Summary 

[35]  In  this  study,  upper  ocean  processes  in  the  tropical 
Pacific  are  examined  in  coupled  models  as  a  means  of 
increasing  our  understanding  of  the  equatorial  Pacific  cold 
tongue  bias.  Our  analysis  indicates  that  the  error  of  ocean 
heat  advection  contributes  significantly  to  the  cold  SST 
biases  along  the  equator.  Most  CGCMs  produce  overly 
strong  total  (horizontal  plus  vertical)  heat  advection  along 
the  equator  (Figure  6c).  This  is  caused  by  overly  strong 
horizontal  advection  of  heat  (Figure  9b)  due  to  overly  strong 
zonal  currents  (Figure  12a)  and  by  overly  strong  vertical  cold 
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Figure  15.  Linear  regression  of  2°S-2°N  averaged  of  20-year  monthly  mean  data  for  (a)  rx  vs  T050 
(Pa  K  *)  and  (b)  vertical  velocity  (0-50  m)  vs  rx  (x  10  4  ms"1  Pa  l)  from  SODA  ocean  analysis  and 
15  IPCC  AR4  CGCMs.  The  data  are  first  converted  to  10-degree  bins  in  longitude  before  meridional 
averaging.  The  period  for  analysis  is  January  1 980-December  1999. 


advection  (Figure  9a)  due  to  overly  strong  vertical  motion 
and  vertical  temperature  gradient  (Figures  1 2b  and  1 2d).  The 
cooling  bias  from  the  total  heat  advection  has  a  magnitude 
similar  to  that  of  surface  heat  flux  (Figure  6b).  Thus,  errors  in 
both  net  surface  heat  flux  and  total  heat  advection  contribute 
to  cold  SST  bias  along  the  equator.  Our  results  suggest  that 
the  biases  in  the  oceanic  processes  discussed  here  contribute 
to  the  eastern  Pacific  cold  tongue  bias. 

[36]  Bjerknes  feedback  analysis  shows  that  almost  all 
models  underestimate  the  strength  of  Bjerknes  feedback  in 
the  central  and  eastern  Pacific  between  160°E  and  110°W 
(Figure  13),  which  may  be  related  to  insufficient  surface 
zonal  wind  response  to  SST  in  the  models  and  an  erroneous 
subsurface  temperature  structure.  An  interesting  result  from 
the  analysis  of  observations  is  that  the  horizontal  heat 
advection,  which  is  dominated  by  its  zonal  component, 
contributes  more  than  the  vertical  advection  does  to  the  total 
advection.  Therefore,  the  total  heat  advection,  instead  of  the 
vertical  heat  advection,  should  be  considered  when  discuss¬ 
ing  about  Bjerknes  feedback. 

[37]  The  overly  weak  surface  wind  response  to  SST  and  the 
resultant  overly  weak  Bjerknes  feedback  in  the  coupled 
models  are  surprising,  since  Lin  [2007]  found  an  excessive 


surface  wind  response  to  SST  for  the  atmospheric  compo¬ 
nent  of  these  coupled  models  forced  by  observed  SST.  Our 
finding  suggests  that  coupling  to  an  ocean  model  signifi¬ 
cantly  reduces  the  excessive  response  of  zonal  winds  to  SST 
anomaly,  which  may  be  caused  by  the  overly  cold  equilib¬ 
rium  SST  along  the  equator  and  the  erroneous  ocean 
subsurface  temperature  structure.  We  hypothesize  that  the 
equatorial  Pacific  cold  tongue  bias  in  the  coupled  models  is 
developed  in  two  stages. 

[38]  Stage  1 :  In  the  initial  state,  both  SST  and  the  ocean 
temperature  profile  are  similar  to  the  observations  in 
CGCMs.  Next  there  is  an  overly  strong  Bjerknes  feedback 
caused  by  the  AGCMs’  excessive  surface  zonal  wind 
response  to  SST  [Lin,  2007],  which  may  be  enhanced  in 
some  models  by  the  excessive  ocean  upwelling  response  to 
trade  winds  (Figure  15b).  Then  the  cold  tongue  becomes 
colder  than  the  cold  tongue  in  the  observations. 

[39]  Stage  2:  Eventually,  when  the  cold  tongue  gets  too 
cold  (Figure  1),  and  the  subsurface  temperature  structure 
becomes  erroneous  (Figure  14),  the  surface  zonal  wind 
response  to  SST  becomes  weaker  (Figure  15a)  and  the 
Bjerknes  feedback  is  significantly  weakened.  This  bias  is 
further  enhanced  by  an  erroneous  ocean  temperature  structure 
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in  most  models.  The  Bjerknes  feedback  then  becomes  weaker 
than  in  the  observations,  and  it  no  longer  contributes  to  the 
excessive  cold  tongue. 

[40]  We  plan  to  test  these  hypotheses  in  future  studies 
using  coupled  model  experiments  initialized  from  the  observed 
SST  and  ocean  structure.  Zhang  and  Wang  [2006]  did  a 
pioneering  study  along  this  line  with  an  emphasis  on  the  role 
of  the  atmospheric  model.  More  work  is  needed  on  the  role 
of  the  ocean  model  and  ocean-atmosphere  feedback  in  the 
growth  of  the  doublc-ITCZ  pattern. 
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